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In their first-principles calculations of the electronic band structure of graphene under uniaxial strain, Gui, Li, 
and Zhong [Phys. Rev. B 78, 075435 (2008)] have found opening of band gaps at the Fermi level. This finding is 
in conflict with the tight-binding description of graphene which is closed gap for small strains. In this Comment, 
we present first-principles calculations which refute the claim that strain opens band gaps in graphene. 



PACS numbers: 73.22.-f, 73.61.Wp, 72.80.Rj 

Gui et ah have used first-principles calculations to inves- 
tigate the effect of planar strain on the electronic band struc- 
ture of graphene, and have found opening of band gaps at the 
Fermi level resulting from arbitrarily small uniaxial strains, 
applied parallel or perpendicular to the C-C bonds. However, 
tight-binding (TB) model on the honeycomb lattice with dif- 
ferent nearest-neighbor hoppings in the three directions has 
been rigorously shown to be closed gap as long as the hop- 
pings satisfy the triangle inequality. 2 ' 3 The closed-gap TB 
model, which contains zero modes, has been further devel- 
oped recently to include the effects of magnetic fields^ and 
corrugations in graphene. 6 The discrepancy between the first- 
principles calculations of Ref. 1 and TB model has already 
been discussed 7 but the reason has not been identified. One 
suggested explanation 7 is that this band gap opening is an ar- 
tifact of density-functional theory (DFT) calculations. How- 
ever, the possibility that the nearest-neighbor TB model is an 
incomplete description has not been ruled outi^ 

We remind that the DFT methods essentially solve single- 
particle Schrodinger equations (Kohn-Sham equations) for ef- 
fective potentials based on the underlying lattice, and the TB 
model solves the same problem in a simplified approximation. 
Therefore, it seems unlikely that qualitative differences exist 
between DFT and TB band structures. We also see indications 
of possible error in Ref. 1. First, Figs. 3(c) and 5(c) of Ref. 1 
show a peculiar peak whose underlying cause is not explained. 
Second, an energy gap is incorrectly ascribed to the TB band 
structure which is then plotted in Fig. 4 of Ref. 1 with large 
symbols that hide the important band crossing. 

In this Comment, we check directly the first-principles 
calculations of Ref. 1 by one of the available DFT codes. 
We used the QUANTUM-ESPRESSO 9 package based on the 
pseudopotential plane-wave method. We obtained the pseu- 
dopotential C.pw91-van_ak.UPF also from Ref. 9 and used 
a kinetic-energy cutoff of 40 Ry, a Monkhorst-Pack fc-point 
mesh of 21 x 21 x 1, and a vacuum separation of 20.5 A along 
the c axis. We chose these parameters as close as possible to 
those of Ref. 1 for a more meaningful comparison^ 

First, we determined the equilibrium lattice constant of 
graphene in the absence of strain. We found a value of a = 
2.464 A, defined in Fig. 03a), which is not significantly differ- 
ent from the 2.4669 A found in Ref. 1 . We then made calcula- 
tions on graphene under uniaxial strain for two special cases, 



for which Ref. 1 has found maximum values for band-gap 
openings. These two cases are (i) 12.2% strain applied par- 
allel to the C-C bond, i.e, along b' direction in Fig. Ufa), and 
(ii) 7.3% strain applied perpendicular to the C-C bond. For 
comparison, recent experimental studies of strain in graphene 
have applied strains of up to 1.3% by stretching or bend- 
ing a flexible substrate, on which graphene was deposited, 
and have measured them by Raman spectroscopy. 11 ! 12 ' 13 Rel- 
atively large strains used in our calculations are more con- 
venient to demonstrate the effects but the conclusions apply 
equally to smaller strains. The uniaxial strain deforms the tri- 
angular lattice of graphene into centered rectangular lattice, 
shown in Fig.Q] Thus for case (i), b' is fixed at 12.2% larger 
value than its unstrained value of ay/3, and the value of a is 
varied until the stress in the a direction becomes vanishingly 
small. Of course, for each choice of a, the positions of the 
atoms must be relaxed until interatomic forces become suffi- 
ciently small. We found Poisson's ratio to be ~ 0.10 for case 
(i). A similar procedure is used for case (ii), with fixing a at 
a value of 7.3% larger than the original value, and then opti- 
mizing b' . Here Poisson's ratio was found ~ 0.14. 

The band structure we obtained for case (i) along the k- 
point path of Fig.[Tfb) is shown in the top panel of Fig. [2] with 
its important portion magnified in the bottom panel. In the 
band plots, we used a regular fc-point mesh of 60 points for 
the entire path and refined it by the addition of 20 extra points 
as shown in the magnified part. The existence of a contact is 
clearly seen between the conduction and valence bands near 
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FIG. 1: (a) (Color online) Graphene honeycomb lattice. It is formed 
from a triangular Bravais lattice with a two-atom basis consisting of 
A\ and B\. Under uniaxial strain the Bravais lattice becomes cen- 
tered rectangular, (b) First Brillouin zone. TKM and TKRS enclose 
the irreducible wedges corresponding to the triangular and centered 
rectangular lattices, respectively. 
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FIG. 2: (Color online) Band structure of graphene under 12.2% strain 
applied parallel to C-C bond (top). Magnified portion of band struc- 
ture near the contact point (bottom). 

K on the YK line of the Brillouin zone. The displacement of 
the Dirac cone along YK toward Y is in agreement with the 
TB description of Ref. 2. Our Fig. 2 is to be compared with 
Fig. 3 of Ref. 1, where they give a value of 0.486 eV for the 
band gap. In our calculation, we find a value of 0.498 eV for 
the energy splittings at K and R in the top panel of Fig. [2] 
which is probably what is taken, in this case, as the band gap 
by Ref. 1, having missed the nearby band crossing. 

The band structure corresponding to case (ii) is shown in 
Fig- El an d must be compared with Fig. 5 of Ref. 1. The main 
difference with case (i) is that here the band crossing occurs 
near the R point on the RS line. This is equivalent to a shift 
of the Dirac point along the YK line away from Y, i.e., in 
the opposite direction to that of case (i). In Ref. 1 a value of 
0.170 eV is given for the band gap for this case. We found a 
value of 0.178 eV for the energy splittings at R and K, which 
is, as in the other case, close to the band gap given in Ref. 1. 



In conclusion, our first-principles calculations establish that 
graphene under uniaxial strain is gapless in agreement with 
the tight-binding model. Our numerical values and the general 
shape of band structures are quite similar to those found in 
Ref. 1 . However, accidental degeneracies in the band structure 
have been disregarded in Ref. 1, and this has resulted in the 
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FIG. 3: (Color online) Band structure of graphene under 7.3% strain 
applied perpendicular to C-C bond (top). Magnified portion of band 
structure near the contact point (bottom). 



appearance of spurious maxima in band gaps as a function of 
strain. 
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